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1. Introduction 


The problem of extending the thermodynamic formalism to systems out of equilibrium is 
a central challenge in statistical physics. Signihcant progress on this topic has been made 
in the past two decades, yielding several remarkably general results such as fluctuation 
theorems This has resulted in a large body of work that extends the notion of 

statistical ensembles to dynamical trajectories [6 -13 , and is closely connected to the 


mathematical theory of large deviations, which deals with probabilities of rare events 
(for a review, see e.g. [^). In particular, one is interested in calculating the statistics 
of time-integrated observables At[x\, where x is a (phase-space) trajectory of duration 

r. 


P[At = a) = V[x]P[x]5{At[x] - a). 


( 1 ) 


Here P[x] is a probability density functional of a path x, and the integration is over 
all possible paths. Examples of A^ include the action functional [^, average drift of 
a Brownian particle 16 , transition history ^ and time-integrated particle current in 


driven diffusive systems, e.g. in exclusion processes 17 


Often, it is of interest not only to calculate the probability density P{At = a) of a 
rare event, but to understand how the particular fluctuation leading to the rare event 
occurred. This amounts to selecting trajectories x that have fixed At[x\ = a, 

P[x, At = a] 


P[x\At = a] = 


( 2 ) 


P{AT = a) ’ 

which resembles the microcanonical ensemble. Unless we are in the low-noise regime, 
there will be many trajectories contributing to the event At = a. Selecting trajectories 
leading to At = a from the original, unconstrained dynamics may prove difficult if 
At = a is a rare event (a large deviation), in which case the unconstrained average (^r) 
is different from a for large T, where (...) is taken with respect to P[x]. Even worse, 
we do not generally expect to find the constrained process to be Markovian, which can 
pose serious difficulties for its analysis. A resolution to this problem in the spirit of 
equilibrium statistical physics, that is attracting much interest of late, is to look at the 
canonical path ensemble instead. 


PJx = 


(esAr^ 


( 3 ) 


This path ensemble is known under several names, depending on the field of study. 
In physics, it has been called s ensemble (related to the choice of letter s for the 
tilting parameter), but also driven or biased ensemble; in rare-event simulations it is 
called (exponentially) tilted ensemble, or Esscher transform of P. It has been used to 
probe rare trajectories in systems with metastable states, namely to study the glass 
transition 18 -25 and more recently to address the problem of protein folding 26 -28 


Recently, Chetrite and Touchette have rigorously proved several remarkable 


properties of the canonical path ensemble 12,13 . First, the canonical path ensemble 
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can be realized by a Markov process (referred to as driven dynamics) in the long-time 
limit. Second, the microcanonical and the canonical path ensembles are asymptotically 
eqnivalent in the limit T —)■ oo in the sense that 


ii„, iioghidn^ 

T^ocT ^ PJx] 


= 0 


( 4 ) 


where the limit is approached almost snrely and P[x\Arp = a] and Ps[x] are given by 
(|^ and ([^ respectively. This important resnlt allows one to study (generally unknown) 
dynamics of constrained systems by studying dynamics of the driven process, which can 
be then analyzed using standard Monte Carlo techniques. The equivalence holds under 


the following three conditions 13,29 


Condition A: P{At = a) must satisfy a large deviation principle with a rate 
function 1(a), 


where 


F(Ar = a) X T ^ oo, 

denotes asymptotic behaviour in the sense that 
1(a) = — lim 4 log P(At = a). 


T—>-oo 


( 5 ) 


( 6 ) 


• Condition B: the large deviation At = a must arise entirely from the interior of 
the interval [0,T], e.g. excluding large deviations at the boundaries 0 or T, 

• Condition C: the rate function 1(a) must be convex at a. 

We note that the rate function 1(a) does not need to be differentiable at a for the 
equivalence to hold. However, if it is then the tilting parameter is unique and is given 
by s = I'(a). Conditions A, B, C are found to be sufficient for the ensemble equivalence 
to hold; whether all conditions are necessary remains to be clarihed. 

In this work, we are interested in cases where the microcanonical and canonical 
path ensembles are not equivalent. The inequivalence of the nonequilibrium path 
ensembles has been recently reported for large current fluctations in the zero-range 
The exactly solvable one-site model revealed that the problem lies in 


30 31 


process 

the inhnite state space. For example, the spectrum of the driven process may become 
gapless or the process may have non-normalizable eigenvectors with respect to the initial 
measure. Beyond this example, the necessary properties of the spectrum required for 
the equivalence of nonequilibrium ensembles is an open problem. 

In equilibrium, the inequivalence of microcanonical and canonical ensembles has 
been extensively studied for various models (for a review, see [^), typically for systems 
with long-range interactions, such as gravitational systems, dipolar or unscreened 
Coulomb systems. The total energy in these systems is no longer extensive in the energies 
of the subsystems, which may lead to the microcanonical entropy that is a nonconcave 
function of the energy. In that case the canonical free energy is not differentiable at some 
point and thus there is no temperature to establish the connection with the canonical 
ensemble . Consequently, these systems show interesting phenomena such as negative 
(microcanonical) specihc heat, slow relaxation, quasi-steady states, etc. 
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This motivates us to pose the following questions. What are the typical situations, 
analogous to e.g. long-range interactions in equilibrium, where the equivalence of path 
ensembles is not expected to hold? Moreover, are there common phenomena arising 
from the inequivalence of nonequilibrium path ensembles? 

which are 


To address these questions we focus here on stochastic bridges 34 36 


obtained by conditioning Markov processes on fixed xq and xt- Thus the observable At 
is given by 

x{T) — a;(0) 


At — 


T 


(7) 


Crucially, stochastic bridges are themselves Markov processes for any finite time T (a 
property that is unlikely to hold for conditioning on more general At) which allows 
us to study their dynamics analytically. We present the solutions of several stochastic 
bridges for which the original, unconstrained dynamics violates one of the conditions 
A or B for ensemble equivalence mentioned above. These exact solutions shed light 
on why ensemble equivalence breaks down. For the examples we consider, ensemble 
inequivalence is typically manifested by condensation-like phenomena: the conditioning 
At = a causes the original process to ‘condense’ its large deviation (to meet the 
conditioning) in a vanishing fraction of the interval [0,T], rather than throughout the 
whole interval as implied by the driven process associated to the canonical ensemble ([^. 

Related condensation phenomena are well known in the problem of sums of 
independent and identically distributed random variables. If the distribution of the 
random variables is heavy tailed then a large deviation of the sum typically occurs 


through a single random variable realising the large deviation of the sum 37 40 


Moreover, recent work has shown that a heavy-tailed distribution is not necessary for 
condensation to occur when there is a further constraint on the random variables in 

This type of condensation is manifested in real-space 


addition to their sum 41,42 


condensation in spatially extended stochastic mass transport models such as the zero- 
range process wherein a single site captures a finite fraction of the total mass in the 
system 43 . In this work we will examine the connection between temporal condensation 


exhibited in discrete-time stochastic bridges and condensation within collections of 
discrete random variables. 

We start by reviewing known results for diffusion bridges previously presented 
in [^: the equivalence of ensembles holds for Brownian motion, but fails for the 
Ornstein-Uhlenbeck process where the large deviation is condensed into a boundary 
effect. We extend this result for temporal condensation to other diffusion processes in 
which there is a form of potential that causes the conditioning xt = clT to be met in a 
non-extensive manner, to be explained in detail later. To stress the importance of the 
condition B for the equivalence to hold, we present a diffusion bridge for a multiplicative 
noise process (the CIR bridge) in which P{At = a) has an exponential tail, but where 
the large deviation At = a is again a boundary effect. As a borderline case, we also 
study the Brownian motion with dry friction 44,45 , for which the path ensembles are 
equivalent, but the connection between the two ensembles is not unique - a case known 
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as partial equivalence 33 46 


After studying diffusion bridges, we look at several discrete-time Markov chains 
driven by heavy-tailed (and thus non-Gaussian) noise. There the rate function is 
formally zero and thus violates the condition A, which can be also related to the 
phenomenon of condensation in stationary states of mass-transport models (for a review, 
see [^). At the end, we use the equivalence (|^ to revisit the interaction-driven 
condensation 47 , and show that it is equivalent to the conditioning of random walk 


trajectories on a large deviation of their local time, which counts the number of returns 
to the origin. Our results should serve as guiding principle for other, more complex 
stochastic systems, for which conditions A or B are generally difficult to examine. 

The paper is organised as follows. We study stochastic bridges for diffusion 
processes in Section and for discrete-time Markov chains driven by heavy-tailed noise 
in Section!^ Discussion and conclusions based on these examples on when not to expect 
the equivalence to hold are presented in Section 


2. Diffusion bridges 

In this paper, we consider a diffusion process dehned by the following (ltd) stochastic 
differential equation (SDE), 


(JjX 

— = b{x,t) + a{x,t)ri{x,t), a:(0) = 0 


( 8 ) 


where rjit) is h-correlated noise, 

{v{t)) = 0, {r]{t)r]{t')) = 6{t-t'). (9) 

A diffusion bridge is a process obtained by conditioning x{t) in (|^ to have a hxed value 
of K at time T[tJ It is an ideal candidate for studying nonequilibrium path ensembles, 
for two reasons. First, the observable At ([^ takes the simple form 

( 10 ) 


. a: aT 
= ^ = ^ = a, 


where we have set K = aT to impose a large deviation. Second, the conditioned process 
itself can be conveniently described by a stochastic differential equation that is similar 
to (j^, but with a modified drift term [I^[M[j3^ . To see this, let p(a:, t|0,0) solve the 
Fokker-Planck equation corresponding to the stochastic differential equation (|^, subject 
to the initial condition p(x, 0|0, 0) = h(x). 


dp{x,t\0,0) d 1 ^ ^ ^ 

— Ft — = “to w [<^ (^.*)p] ■ 

The probability density g{x,t) for the stochastic bridge is given by 

p(A:,T|x,f)p(x,t|0,0) 


g{xA) = 


p{K,T\0,0) 


( 11 ) 


( 12 ) 


f This value is traditionally set to 0; here we adopt the name “bridge” for a diffusion process on [0, T] 
conditioned on x{T) = K, irrespective of the value of K. 
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Here, p{K,T\x,t) solves the backward Kolmogorov equation 

^^p{K,T\x,t) = ( 13 ) 

subject to the final condition p{K^ T\x, T) = 6{K — x). 

Now it is straightforward to obtain 

dg{x,t) _ p{K,T\x,t) dp{x,t\Q,Q) p{x,t\Q,Q) dp{K,T\x,t) 

dt p(K:,T|0,0) dt p{K,T\0,0) dt ' ^ ^ 

Then substituting ( IIpK ) and regrouping terms yields, after some calculation, the 
following Fokker-Planck equation for g{x,t) 


dg{x,t) d 


1 


dt 




2 dx"^ 


where CT{x,t) is given by 


d 


crix^t) = b{x,t) + a (x,f)—lnp(iF,T|a;,f), 

dx 


(15) 


(16) 


and the subscript T emphasises that CT{x,t) depends explicitly on T. Using (15), we 
arrive at the following stochastic differential equation for the diffusion bridge which we 
denote here and in the following by y{t) 


^ ^CT(,y,t) + a(y,t)7](y,t), 


( 17 ) 


where the expression for CT{x,t) is given in (16). Note that the equation 0 is an 
exact equation for the conditioned process which is a diffusion with (non-homogeneous) 
drift given by (16) and noise width a. This fact allows us to study the microcanonical 
dynamics, provided we can solve the backward Kolmogorov equation for p{K,T-, x,t). 
If so, the equation (17) can be then solved numerically using e.g. an Euler-Mayurama 
scheme. Moreover, if p{K, T]x,t) has a simple enough form we can compute exactly 
various quantities for the conditioned process. Next, we will present several examples 
of diffusion bridges, for which there is an explicit expression for g{x,t) and thus for 
CT{x,t). 


2.1. Brownian bridge 


For the Brownian bridge discussed e.g. in Ref. 13 , we have 

b{x,t) = const. = /i, a{x,t) = const. = a, (18) 

which also includes a special case of p = 0 called the Wiener bridge. The solution to 
the corresponding Fokker-Planck equation subject to p{x, 0) = 6{x) is given by 

1 


p(a;,t|0,0) = 


-_e 2ta^ , 


(19) 


a\/2TTt 

Because the process is both space- and time-homogeneous, the expression for p(iF, T|a;, t) 
is the same as for p{K — x,T — t\0,0), yielding 

1 _ [K-x-fi{T-t)]^ 

_Q 2(r-po-2 


p{K,T\x,t) = 


a 


y2ir(T - () 


( 20 ) 
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For K = at, the stochastic differential equation for the Brownian bridge is thus given 
by 

dy aT-y 

* = pi) 

where we denoted the bridge by y{t) to distinguish it from the unconstrained process. 
We note that the drift (aT — y)/(T — t) of the bridge process is actually time and space 
dependent, but as we shall see, for large y becomes independent of T and ~ a. Also 
note that equation (21) is linear in y and therefore its solution can be written explicitly, 

dW, 


y{t) = at + a{T — t) 


'0 


T- s’ 


( 22 ) 


where Wg is Wiener process. The process y{t) in (22) is Gaussian, whose mean and 
variance can be easily evaluated from 


and read 
{y{t)) = at, {y{tf) - (2/(t))^ = 


1 - - 
T 


(23) 



Figure 1. Brownian bridge starting at ?/(0) = 0 and ending at y{T) = aT for /r = 0, 
CT = 1, a = 10 and T = 1, calculated at discrete time intervals of size e = 0.001; solid 


lines are {y) and {y) ± 3[(?/^) — where (y) and (y^) — (y)^ are given in (23). 


A sample trajectory for the constrained process, obtained at discrete time intervals of 
size e = 0.001, is presented in Figure 

The equivalence between the Brownian bridge and the corresponding driven process 
has been established in Ref. Essentially, one looks at the probability density 


F ( At = ^ = a I = 


e"^ 2,72 


P{At = a) = P{x{T)/T = a) for the unconstrained dynamics, which from (19) reads 

The corresponding rate function is given by 

- i-^r 


(24) 


/(a) = 


2 ( 7 = 


(26) 


and is differentiable, yielding s = P{a) = as the tilting parameter. One then 


constructs ‘driven dynamics’ using the generalised Doob transform 13 , to obtain that 














the drift of the driven process in the large T limit is exactly a, as for the Brownian 
bridge. 

The Brownian bridge is an example where the (global) conditioning modihes 
trajectories locally in the interior of the interval [0,T]. The next example - that of 
the Ornstein-Uhlenbeck bridge - is rather different, in the sense that the large deviation 
that fulhls the conditioning is concentrated at the end of the interval. 


2.2. Ornstein-Uhlenbeck bridge 


Next we review results for the Ornstein-Uhlenbeck process, also studied in Ref. 13 , for 
which 

b{x,t) = — x), a{x,t) = const. = a, (26) 

where 6 > 0. The solution to the corresponding Fokker-Planck equation subject to 
p{x, 0) = 6{x) is given by 

[x-/j(l-e-'^*)]" 


p(a;,f|0,0) = 


■2cr2(l-e-2»‘)/(2e) 


Similarly, the expression for p{K,T\x,t) is given by 

^K—^—(x— 

p(K,T\x,t) = 


(27) 


2t72[i_e-20(T-t)]/(26)) 


a 


y2ir|l-e-2‘'(’’-<)|/(2»)' 


(28) 


For K = at, the stochastic differential equation for the Ornstein-Uhlenbeck bridge 13 
denoted y{t), is thus given, from (|I6|) and (17), by 


dy cosh(6*(T — t)) 6{aT — p) 

smh(«(T - ()) + smh(«(T - t)) + 

This equation is again linear in y and its solution can be written explicitly as 

sinh(6'(T — t)) 


(29) 


, . , ^ , sinh(6'f) 

y{t) = {aT - p) . ’ + p 

smh(fc^i ) 

rt 

-\- asmh{9{T — t)) 


1 - 


sinh (6'T) 


+ 


dhU, 


/o sinh(6'(T — s)) ’ 

As a Gaussian process, y{t) is fully determined by its mean and variance which read 

sinh(6'(T — t)) 


(30) 


, / ^ N sinh(6'f) 


1 - 


sinh(6'T) 




ivW ) - (»(<)> = 


(31) 

(32) 


sinh (6*T) 

A sample trajectory for the constrained process, obtained at discrete time intervals of 
size e = 0.001, is presented in Figure]^ 

For 6t, 9T S> 1, the average of y{t) can be approximated by a simpler expression 
{yit)) ^ {aT - + p{l-e-^y , 0f,0T>l. (33) 
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Figure 2. A sample trajectory of the Ornstein-Uhlenbeck bridge starting at j/(0) = 0 
and ending at y{T) = aT for 0 = 10, ^ = 0, ct = 1, a = 10 and T = 1, calculated at 
discrete time intervals of size e = 0.001; solid lines are (y) and (y) ± 3[{y^) — 
where (y) and (y^) — (y)^ are given in (311. Vertical line is placed at T — 1/9 and 
denotes a characteristic time scale over which the large deviation occurs. 


In contrast to the Brownian bridge y{t) in Figurej^where {y{t)) = at, here {y{t)) has the 
same value ~ /r as in the unconstrained process, except for t approximately 1/9 away 
from 0 and T. Put differently, the conditioning to reach aT at time T has no effect on 
the stochastic dynamics in the interior of [0,T], except for a small fraction of time close 
to T, which goes to 0 as T —)■ oo. This boundary effect is not recovered by the driven 
process, which has the same drift as the unconstrained process [^. The driven process 
correctly describes the conditioned process in the interior of [0,T], but not at the right 
boundary. 

Here we give an intuitive explanation for this effect. The time-integrated observable 
Arp, such as At = [a;(T)—a;(0)]/T, is nothing more than a sum (or an integral) of random 
variables. Loosely speaking, the large deviation At = a is an interior effect if all of these 
random variables can be changed locally to achieve the large deviation. To see that this 
is not the case in the Ornstein-Uhlenbeck process, we consider its discrete-time Markov 
chain analogue, called the autoregressive process of order 1, or AR(1), 

Xt = Xt-i + 9(y — Xt_i) + arjt) = 0 (34) 


where we assume that 0 < 6* < 1, and rjt are independent and identically distributed 
Gaussian random variables with the mean {rj) = 0 and variance (r^^) = 1. The stochastic 
recurrence equation (34) can be iterated yielding 

T 


Xt = y[l - (1 - 9f] + ^^^(1 - (35) 

t'=i 

In principle, all yp contribute to a large deviation of their sum. However, the weighting 
factor (1 — 9)’^~*'' will make only few of them contribute to the sum, ones that are close 
to T. In this sense, we may say that the boundary effect is due to the fact that Xt 
is not extensive in T - increasing T will not increase the number of random variables 
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rjt! that contribute substantially to Xt- 
following expression, 


From (27) we also get for P{x{T)/T 


a) the 


_ ea^P 

P = a) oc e , T —)■ oo, (36) 

so that the rate function I (a) dehned with respect to the limit T —)■ cx) in ([^ is formally 
inhnite. However, it is important to emphasise that the main reason the path equivalence 
(|^ does not hold here is really the boundary effect described above. In the next example, 
we will study a process for which P{x(T)/T = a) has an exponential tail, but the 
equivalence does not hold because of a similar boundary effect. 


2.3. Cox-Ingersoil-Ross (CIR) bridge 

Here we study a process that has the same drift as the Ornstein-Uhlenbeck process, but 
with a different, state-dependent diffusion coefficient 

b{x,t) = 6{ij, — x), a(x,t) = a\/x, x>0. (37) 


This process is known as the Cox-Ingersoll-Ross model or the CIR process, and belongs 
to a class of square-root diffusions. The CIR process is a popular model for the evolution 
of interest rates in finance 48 . Remarkably, the propagator p{x,t\x',t') of the CIR 


process is known in a closed form 49 


, I , 26 / X \ 

A6^/ xx' 


.q/2 29(3; + 3:'e-^(*-*') 

' g a2[l-e-9(‘-‘')] 


X In 


(38) 


cT^[l — e *')] 

where q = {29p/a‘^ — 1) and Iq{x) is modified Bessel function of the first kind. The 
expression for p{aT,T\x,t) can be then used to calculate the drift ct^xR), and the 
corresponding SDE for the bridge can be then integrated numerically. However, we note 
that since x must be non-negative at all times, the standard Euler-Mayurama scheme 
is not appropriate for numerical integration. Here, we used modified Euler scheme 
for SDEs with square-root diffusion coefficient described in Ref. 50 , which has strong 
convergence. A sample trajectory for the constrained process, obtained at discrete time 
intervals of size e = 0.001, is presented in Figure Sample trajectories of the CIR 
process look similar to those for the Ornstein-Uhlenbeck process, the difference being 
that larger values of x are more stochastic, which is due to the factor -y/x multiplying 
white noise term. On the other hand, the limit 


1(a) = lim —InP 
^ ^ T-s-oo T 


T 


= a \ = 


200 

7^ 


(39) 


is finite, which emphasises the importance of the condition B for establishing the 
equivalence even when the condition A is satisfied. 

The examples of sections 2.2 and 2.3 illustrate the phenomenon of ‘temporal 
condensation’ (the name has been suggested in (^), whereby a large deviation is 
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t 


Figure 3. A sample trajectory of the Cox-Ingersoll-Ross bridge starting at 2/(0) = 1 
and ending at y{T) = aT for 0 = 10, /i = 1, ct = 1, a = 10 and T = 1, calculated at 
discrete time intervals of size e = 0.001. 


localised to a small fraction of the time interval T (that vanishes in the limit T —)■ cx)), 
rather than distributed over the whole time interval as for the Brownian motion in 
Section 12.11 

Next, we consider a borderline case - a Brownian motion with dry friction - which 
interpolates exactly between the Brownian bridge and the last two cases, yielding a finite 
fraction of the time interval T that is spent around the mean value of the unconstrained 
process. 


2.4- Brownian bridge with dry friction 


Brownian motion with dry (or Coulomb) friction is governed by the following SDE: 


m— = —yn — A/ 7 Sgn(n) + maT]{t), 

(JjL 


(40) 


where 7 is a friction coefficient and A/? is a threshold force. Here we consider only the 
inviscid case 7 = 0, which was first proposed by de Gennes to describe a particle on a 


vibrating plate 44 ; for 7 7 ^ 0, see work by Touchette, Van der Straeten and Just 45 


By introducing A = Ap/m and renaming x = v, we can write (40) with 7 = 0 as 

(41) 


dx A / N 

— = -Asgn(a;) 


ar], 


which identifies h{x,t) = — Asgn(x) and a{x,t) = const. = a. The solution p{x,t\x',0) 
to the corresponding Fokker-Planck equation is known 45 and is given by 

g-A2t/(2<T2) 

p{x,t\x', 0 ) = 


g-A(|a:|-|a;'|)/a2g-b^ 


A- 


j-2A|a;|/o-2 

2^2 


■\/2'KaH 




V2aH 


(42) 


where erf(a;) is the error function. Since the process is time-homogeneous, we can write 
p{x, t\x', t') = p{x,t — t'\x4 0 ), yielding the following expression for the drift ct{x, t) from 
















( 16 ), 
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aT—x 


aT(x-\-\x\) 


Asgnfa:) + 

CT{x,t) = — Asgn(x) H- 


Asgn(a;)e 


(T-t), 


aT(x+\x\) 


(43) 


HT-t) e'^Terfc(JT) 


1 ^ V 

where erfc(a;) = 1 — erf(a;) and Jt is given by 
jf x\ {a - A)T + \x\ + At 

^ ^ V2a2(T - 1) ^ ^ 

We can easily check that for A = 0 the drift term CT{x,t) rednces to (aT — x)/{T — t) 
which is the case of the Brownian bridge (21). 

First let us consider ct{x, t) for small t in the limit when T —)■ cx) for a; ~ 0 , in which 
case Jt —)• — C )0 for a < A and Jt —)■ C )0 for a > A. Using the asymptotic expansion of 
the complementary error function, erfc(a:) ~ [sgn(x) — 1] + exp{—x^)/{^/nx) for large 
|a:|, we obtain 

-Asgn(a:), a < A 

crix^t) = { aT-x ■ (45) 


T-t ’ 


a > A 


This behaviour of the drift CT{x,t) implies that, in the beginning, the process starting 
at x(0) = 0 behaves in the case a < A as the original, unconstrained process which stays 
close to X = 0, or in the case a > A as a Brownian bridge which on average ascends 
with drift a. 

Let us now consider the case a < A in which case the process begins like the 
unconstrained dry friction process. Consider the limit when both T, t —)■ cx) with r = t/T 
for X ~ 0. Then Jt —)■ —oo for r < 1 — a/A and Jt —)■ C )0 for r > 1 — a/A. As before, 
we obtain asymptotically 


CT{x,t) 


—Asgn(x), 
aT — X 
T-t ’ 


r < 1 — a/A 
r > 1 — a/A 


(46) 


This behaviour implies that for a fraction of the duration r = 1 — a/A the process 
behaves as the original, unconstrained process which stays close to x = 0. Then for 
later times t > {1 — a/A)T the process behaves as a Brownian bridge which ascends 
to the end point x = aT. Now in this regime we have x ~ (f — tT)ct which yields 
X ~ A{t — tT) and ct — A. Thus in the Brownian bridge section of the trajectory the 
average drift is A. Note this drift is larger than a. 

Particular realisations of these two cases a < A and a > A are presented in hgures 
l^a) and|^b) respectively. 

Despite the fact that the constrained process for a < A behaves in the same way 
as the unconstrained process for a macroscopic fraction of time T, we can check that 
the large deviation principle (|^ holds. Indeed, from the solution of the Fokker-Planck 
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Figure 4. A sample trajectory of the Brownian bridge with dry friction starting at 
y{0) = 0 and ending at y{T) = aT for (a) a = 0.5 < A = 1 and (b) a = 1.5 > A = 1.0, 
calculated at discrete time intervals of size e = 0.2; other parameters are T = 1000 and 
cr = 1. In (a), red line denotes the expected beginning of the ascent at t = (1 —o/A)T. 


equation in (42) we find that P{At/T = a) satisfies the large deviation principle in (|^ 
with rate function 


I{a) = 


2A|a| 




|a|+A)' 

2(t2 


|a| < A 


|a| > A. 


(47) 


One can also easily show that /(a) is differentiable everywhere except at a = 0, so that 
the tilting parameter s is given by s = I'{a) for a ^ 0. However, because I{a) is linear 
for |a| < A (see Figure]^, the tilting parameter s = 2A/(T^sgn(a) depends only on the 
sign of a and not on its absolute value. In this situation, the driven process with the 
tilting parameter s = 2A/(j^sgn(a) corresponds to a range of constrained processes that 
all have |a| < A: the role of a is to parametrise the fraction of time r (1 — a/A)T 
spent around the steady state of the unconstrained process before it typically starts 
to ascend with the drift m This can be thought of as phase coexistence between a 
drift-free process and a process with drift A, both lasting hnite fractions of the duration 
T. This phenomenon is known as partial equivalency the hallmark of which is exactly 


this kind of phase coexistence 33,46 


Let us conclude the analysis of diffusion bridges in this Section by mentioning that 
the last three diffusion processes - the Ornstein-Uhlenbeck process, the CIR process 
and the Brownian motion with dry friction - are somewhat special, in the sense that 
these processes respect detailed balance. This fact can be then exploited to understand 
why the trajectories depicted in Figures]^ and ^a) behave initially as if they are not 
affected by the conditioning = aT. 

To this end, let us consider a diffusion process x{t) on the interval [0, T] that starts 
at x(0) = aT, and has a time-independent drift b{x) and a time-independent diffusion 

f Technically speaking, this will be true only in the limit T —>■ oo. For a large, but finite T, the straight 
line in /(a) may have a non-linear correction that selects a particular a. 

f We also note that the rate function I(a) is not differentiable at the point a = 0, which means that for 
0 = 0 the tilting parameter s is not unique. This fact has also been referred to as partial equivalence 


in some works 51 . 
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Figure 5. Rate function I(a) describing large deviations of P{x{T)/T = a) oc 
exp(—r/(a)) for Brownian motion with dry friction; other parameters are A = 10 and 
a = 1. Dashed lines separate linear dependence on a (red lines) where the equivalence 
of nonequilibrium path ensembles is partial, from quadratic dependence (blue lines) 
where the equivalence is full. 


coefficient a(x). Let us assume that p{x,t) is solution of the corresponding Fokker- 
Planck equation with initial condition p{x,0) = 6{x — aT), so that the process x(t) 
starts at aT at t = 0. The detailed balance condition states that 


p{x, T - t\y, 0)ps{y) = p{aT, T - t\x, 0)ps{x), 
where Ps{x) is the corresponding stationary distribution, 

N 


Ps{x) = 


pirAx'h{x')la{x'f 


(48) 


(49) 


aixY 

and M is the normalisation constant. 

Next, we consider the time-reversed process x{t), dehned as x{t) = x{T — t) 
This process is known to be again a diffusion process, governed by the following (Itb 


stochastic differential equation [52 |53 
dx - 

— = b{x,t) -b a{x)p, 
where b{x,t) is given by 

b{x, t) = —b{x) + 


d 


p{x, T — t\aT, 0) dx 


[a‘^p{x, T — t\aT, 0)]. 


(50) 


(51) 


(52) 


Using the detailed balance condition in (48), we can rewrite 6(x, t) in (51) as 

— 0 d(j‘^ d 

b{x, t) = -b{x) + lnp(aT, T - t\x, 0) -b —-b a‘^—\nps{x). 

ox ox Ox 

The last term can be written as 2b{x) — da'^/dx so that the hnal expression for b{x,t) 

is given by 

“ d 

bix, t) = bix) + cr^—lnp(aT, Tlx, f), 
dx 


( 53 ) 
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where we used the fact that p{x,t\x'^t') = p{x,t — t'|x',0). Notice that b{x,t) is the 
same as the drift CT{x,t) derived previously for a diffusion bridge, see (16). We have 
thus shown that the process conditioned to reach aT at time t = T is in fact the time- 
reversed process of the one that starts from aT, the only difference being that the end 
point of the time-reversed process x{T) is unconstrained, rather than fixed as for the 
diffusion bridg^ Starting from aT at time f = 0, the Ornstein-Uhlenbeck process and 
the Brownian motion with dry friction both relax to their stationary distributions, which 
explains why the trajectories depicted in Figure and [^a) are initially not affected by 
the conditioning x{T) = aT. 

It is noteworthy that we can extend this result further to simulate any diffusion 
bridge derived from a process that respects detailed balance, but whose full time- 
dependent solution of the Fokker-Planck equation is not known explicitly. The idea, 
proposed in Ref. 54 , is to simulate two processes, one that starts from a;(0) = 0 and 
the other that runs backwards in time and starts from a:(0) = aT, until they intersect at 
some later time f] if 0 < < T, we can construct a diffusion bridge y{t) as a piecewise 

function defined by y{t) = x{t) for t < t' and y{t) = x{t) ioi t' < t < T Finally, let us 
emphasise that the condensation phenomenon studied in this Section is by no means 
due to the detailed balance only. Rather, it is caused by a cost of having a large value 
of Xt for all 0 ^ t ^ T, that is associated with a confining potential R(^) (defined 
such that b{x) = —V'{x)) that grows to inhnity faster than linearly. For example, an 
overdamped Brownian particle in a periodic and bounded potential (with or without 
external driving that breaks detailed balance), conditioned on Xt = aT will show no 
condensation phenomenon 


3. Random walk bridges 

So far we have presented diffusion processes, in which the noise was always Gaussian. 
In this Section we analyse continuous-space random walks in which the noise is not 
necessarily Gaussian. Although such processes typically converge to diffusions, some 
rare events that violate condition A are not captured in this limit. In particular, we are 
interested in a simple random walk on a real line, defined by the following stochastic 
recurrence equation 

Xt = Xt-i + r]t, Xq = 0, (54) 

where r]t, t = 1... ,T are independent and identically distributed random variables with 
a common probability density (p{rit), not necessarily Gaussian; we assume that {rjt) = fi 
and {rj'f) — {rjt)"^ = As before, we are interested in the bridge process, obtained by 
conditioning Xt on fixed value of Xt = aT, which amounts to fixing the value of the 

f That is not so important here, because we used initial conditions that are also stationary points 
(attractors) of the noiseless equation. 
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sum 


XT = '^r]t = aT. 


(55) 


t=i 


The probability density P[X\Xt = aT] of a path Xt conditioned on Xt = aT can be 
written as, 


P[X\Xt = aT] = 


P(Xr = aT) 


P[p(A',|AVi)i5(A'T-fflr) 


P(Xt = aT) y-1 


J]4j(A,-A,_ i)i(Ar-or) 


P (HUpih = aT] 7. 


n »’(’(<)'* “ aT 


( 66 ) 


,t=l 


We note that the final expression for P[X\Xap = aT] has the same factorised form as 
the steady-state probability in the well-studied mass transport models (for a review, see 
(i^), where t denotes coordinate in space, and rjt is the mass at the site t. A distinctive 
feature of these models is the phenomenon of real-space condensation, whereby for a > 
there is a single site that carries a macroscopic fraction (a — fi)T of the total mass aT, 


while the rest of the sites have a typical mass of 0 ( 1 ) with mean /r |^,55 ; for a rigorous 
analysis, see |M|[5^[5^ . This situation occurs only if the underlying steady-state weight 
(p{r]) is heavy-tailed (^[^, in the sense that 


(57) 


dT]ip{ri)e'‘'^ = oo for fc > 0 , 


which is to say that the moment-generating function of rj does not exist. Examples of 
heavy tails are a stretched exponential ^{rj) cx exp(—ar^") with a < 1 or a power law 
(p{ri) cx A/rf with b > 2. In all these cases, P{Xt = aT) takes the following form in the 
limit of large T 


P{Xt/T = a) = T^p{a — fi), T —)■ oc, 


(58) 


SO that the large deviation principle (condition A) does not hold. The particular form of 
P{Xt/T = a) tells us that the event Xt/T = a is realised by a single random variable 
taking a large value; this random variable can be any of the T random variables rjt, 


t = 1,... ,T, hence the prefactor T. An example of the random walk bridge (54) driven 


by a heavy-tailed noise is presented in Figure]^ for (a) a < /i and (b) a > ja, the latter 
showing a distinctive jump of size of (a — /i)T 

Recently, we reported a condensation transition for independent and identically 
distributed random variables rjt, t = 1,... ,T constrained to have fixed values of both 
M = Y^=i Vt = pT and V = = ST 41,42, where p 7 ^ 1 is a parameter. In 

this situation, there is no need for the distribution of to be heavy-tailed - instead 
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Figure 6. Random walk constrained on the fixed value of Xt = 

T = 1024 and (a) a = 0.2 < fj, and (b) a = 1.5 > /r. Here (fi{x) = a/^x + l)"'*'^, 
where a = 5/2 (the Lomax distribution), which has the mean of /r = 1 /(q! — 1) = 2/3. 
Dashed line in (a) has a slope of a = 0.2, compared to dashed lines in (b) which have 
a slope of /r = 2/3; in (b), the size of the jump is (a — /r)T = 5/6r. 


condensation is achieved through the second constraint, which changes the bare (light¬ 
tailed) probability density to a heavy-tailed one. For p < 1, the condensation transition 
happens for 5 > 6c, where 6c is the mean of the effective heavy-tailed distribution, so 
that the value of V falls in the large deviation regime. 

Applied to the random walk in (54), the second constraint for p < 1 takes the 
form of At = it is often called realised power variation, or realised 

quadratic variation specifically for p = 1/2 [^,59 and measures the variation of the 
trajectory. However, as we showed in Ref. [S,42 , the condensation in this context may 
be viewed as a finite-size effect for p < 1 , in the sense that the joint probability density 
behaves as P(M,V) ~ g-Ti{p)+o{T~tp)^ where the correction due to the condensation is 
of sublinear order in L; here 1 < 7 < 1 /p is related to the tail of the probability density 
in the case <p(p) ~ for large 7 . This further means that the equivalence of path 

ensembles is actually restored in the limit T —)■ 00 , in which the size of the jump 
relative to T goes to zero. On the other hand, for p > 1 the jump is of 0{T) and implies 
ensemble inequivalence in the limit T —)■ 00 . 

An example of a random walk constrained on the fixed values of Xt = Vt = o,T 
and V = Ylt=i Vt — is presented in Figure for (a) 6 < 6c and (b) 6 > 6c, for an 
exponentially distributed noise rjt. While the condensate is visible in the noise variables 
Tjt (lower figure in Figure]^), it is only a minor jump (that scales sublinearly with T) 
in the overall sample trajectory Xt. 


4. Conclusions 

The s ensemble approach holds great promise for studying large fluctuations in 
nonequilibrium systems, regardless of whether they obey detailed balance or not. Its 
connection to the corresponding constrained dynamics leading to a large fluctuation is 
often taken for granted, and is expected to hold in majority of cases. The connection has 
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rjt is exponentially distributed for which Sc can be computed analytically and reads 
Sc = = 2. The second plot in (b) shows a condensate at time t = 621 of size 


\/S — ScL. 


been recently established rigorously, based on conditions rooted in the large deviation 
theory. 

In this work, we presented several constrained stochastic systems in which one or 
two of these conditions are not met, and we studied whether they are equivalent to 
their corresponding driven processes in the s ensemble approach. For a one-dimensional 
stochastic variable Xt, we looked at large deviations of its time-integrated speed on the 
interval [0, T], which amounts to conditioning on the hxed value of Xt- Such constrained 
stochastic processes, called stochastic bridges, are particularly convenient for the present 
study because their dynamics is Markovian and can be constructed exactly for the whole 
interval [0, T], and not just in the limit of large T. As a main result, we showed several 
examples in which the constrained and driven dynamics are not equivalent in the limit 
of large T. Notably, this is manifested by condensation-like phenomena, in the sense 
that to meet the conditioning, the constrained process changes only a small portion of 
the dynamics. We have found essentially two types of condensation phenomena in these 
examples, both related to anomalous large deviations. 

The hrst type of condensation is where a large deviation is not realised in the 
interior of the interval [0, T], but is rather a boundary effect realised at time T. One such 
example, the Ornstein-Uhlenbeck bridge, was hrst reported in Ref. [^. In the present 
study we analysed several other diffusion bridges and argued that a similar boundary 
effect leading to ensemble inequivalence is expected whenever a large deviation of Xt 
is ‘penalised’ by a conhning potential. In addition, we presented a borderline case of 
Brownian motion with dry friction, where constrained dynamics is modihed on a hnite 
fraction of the interval that scales linearly with T, which presents an example of partial 
equivalence. The second type of condensation is where a large deviation is realised by 
a single random variable, and is due to heavy-tailed probability distributions; the same 
phenomenon (in space, rather than in time) has been observed in the stationary states 
of mass-transport models such as the zero-range process. 
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The presented examples are by no means an exhanstive list of processes that 
exhibit condensation and ineqnivalence of noneqnilibrium path ensembles. In fact, the 
eqnivalence established in Refs. provides a powerfnl tool for nnderstanding 

other condensation phenomena whose mechanism of condensation is not obvions. Here 


we briefly mention the case of interaction-driven condensation reported in Ref. 47 


in which the original zero-range process was generalised to inclnde hopping rates that 
depend not only on the departnre site, bnt also on its immediate environment, cansing 
correlations between neighbonring random variables. The stationary probability of this 
process still admits a factorised form, which reads 


-P({mi}) cc /nij - pL 


(59) 


2=1 


where the delta fnnction ensnres that the total mass is conserved. In Ref. 
following choice for w{mi,mi+i) was nsed 


47 , the 


w{m, n) = exp [—J|m — n|] ■ exp 




m,0 


( 60 ) 


where 6m,n stands for the Kronecker delta fnnction, and J and U are constants. When 
inserted into (59), the hrst factor in (60) contribntes to the probability of a random 
walk path, which becomes apparent by making the following change in notation: i —)■ t, 
rrii —)■ Xt and w{mi,mi+i) —)■ ip{Xt+i\Xt), where ip is the same as in (56). The extra 
factor in (60), when inserted in (59), appears in the path probability P[X] as a factor 
of exp(t/TL 7 ’), where Lt is given by 


Lt-7^ 


Xtfl- 


(61) 


t=l 


The time-integrated observable Lt measnres the nnmber of returns to the origin and 
is a discrete-time analogue of local time. The hnal expression for the path probability 
P[X] thus has a hard constraint on the sum which is the total area under the 

trajectory, and the tilting factor of exp{UTLT) related to conditioning on local time. 
The latter conditioning ensures that the number of visits to the origin scales with the 
total time T, so that the whole process is a series of random walk excursions (recurrent 
trajectories that stay positive), conditioned on a hxed value of the total area and on 
the total duration T. One can show that these two constraints are responsible for the 
condensation. Details of this calculation, which is similar to the one for the constraint- 
driven condensation 41,42 , will be presented elsewhere. 
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